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Abstract 

We consider the possibility that quark and gluon distributions in the medium 
created in high energy heavy ion collisions may be modified by a power law 
tail at energies much higher than the temperature. We parametrize such a tail 
by Tsallis distributions with an exponent motivated by phenomenology. These 
distributions are characterized by an effective temperature scale that we assume 
to evolve in time like the temperature for thermal distributions. We find that 
including such a tail increases the rates for photon production and significantly 
delays the emission times for photons of a fixed energy. We argue that these 
effects are sufficiently large that they should be able to account for photon yields 
and flow patterns seen in LHC and RHIC experiments. 


1. Introduction 

Theoretical calculations of thermal photon emission from the quark gluon 
plasma [ll[2l|3llll[5l|6l|7l[8l|9l [TOl [HI [12] have difficulties describing the 
10 experimentally observed photon yield and elliptic flow in heavy ion collisions 
da [n [m [Hj. Calculations result in photon transverse momentum spectra 
that are too shallow and approximately a factor of 4 smaller than the ob¬ 
served yield. Further, the photon elliptic flow V2 is computed to be signifi¬ 
cantly lower than the observed V2- The small computed photon V2 is easily 
15 understood because photons are predominantly produced at early times (which 
correspond to high temperatures), when flow has not yet been built up. Vari- 
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2 REVIEW OF PHOTON PRODUCTION IN 1+lD HYDRODYNAMICAL 

MODELS 


ous scenarios have been suggested in the literature to explain the discrepancy 

[lillIHlIiniliniEQlETlEaEslEllIZS]. 


In this paper we study the effect of deviations from thermal quark and 
20 gluon distributions in the plasma. In particular, we show that non-thermal 
tails, implemented by replacing Bose-Einstein and Fermi-Dirac distributions by 
the corresponding Tsallis distribution, have a significant effect on photon yields 
and typical emission times: As opposed to most other mechanisms that have 
been suggested to resolve the thermal photon puzzle, photon emission times are 
25 delayed and at the same time production rates are increased. This could resolve 
the puzzle by generating both increased photon yields and elliptic flow. 

The paper is organized as follows: We show that the photon yield in heavy 
ion collisions is sensitive to the high momentum region of the distribution in 
Section]^ and demonstrate that enhanced tails can lead to both an increase of 
30 the photon yield and a delay of photon emission in Section We then present 
a detailed quantitative calculation of the photon yield from a one-dimensionally 
expanding plasma with Tsallis quark and gluon distributions in Section]^ We 
conclude in Section [5] 

2. Review of Photon Production In 1+lD Hydrodynamical Models 

Photon production in heavy ion collisions is dominated by photons at a 
momentum scale k ^ AT. To see this consider the formula for the rate of 
photon production from a thermalized quark gluon plasma [261127j 



( 1 ) 


35 where E^ is the photon energy, T the temperature of the medium, and Og and 
a the strong and electromagnetic coupling, respectively. 

To compute the photon yield from an expanding ideal gas in I+I dimensions, 
we integrate this rate over time using the time dependence of the temperature. 


which for an ideal gas of relativistic quarks and gluons is given by ^ 


t/to = T^/T^ . 


( 2 ) 
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3 REPLACING THERMAL DISTRIBUTIONS BY TSALLIS 

DISTRIBUTIONS 


The four volume is 

J d^x = At J tdt (3) 

where At is the transverse area of the interaction region. Since 

, dT 1 

tdt(x—7^, ( 4 ) 

upon inserting the rate at a given T into the integral over the expanding system, 
we get 

exp[-E,/T - 41n(r/To)] . (5) 

By determining the stationary point of the exponent we obtain the typical emis¬ 
sion energy for photons to be 

E^ ~ 4T (6) 

up to logarithmic corrections. (The best way to do the stationary phase distri¬ 
bution is in logarithmic coordinates, y = ln(T/To) where dT/T = dy.) 

This emission occurs when we are in the tail of the exponential distribution. 

40 3. Replacing Thermal Distributions by Tsallis Distributions 

Let us now suppose that the quark and gluon distributions were not purely 
Bose Einstein or Fermi Dirac distributions. Since emission occurs at somewhat 
hard momenta, we might expect the distributions would be well approximated 
by an exponential with a power law tail. An example which has this property 
is the Tsallis distribution. 


f{E) = [l + E/{aT)]-\ (7) 

where E is the quark or gluon energy, and a is a free parameter that determines 
the power law at large E. For E aT this distribution is well approximated 
by an exponential, / ~ but for large energy it goes as / ~ {E/aT)~°-. 

For a thermal distribution undergoing 1-1-1 dimensional Bjorken expansion, the 
integral 

P = j d^pf ^ 1/t, ( 8 ) 
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3 REPLACING THERMAL DISTRIBUTIONS BY TSALLIS 

DISTRIBUTIONS 


as it should for a non-interacting gas. For proton-l-proton collisions, the mea¬ 
sured distribution of produced charged particles is approximately a Tsallis dis¬ 
tribution with o « 6. In order that the Tsallis distribution and the thermal dis¬ 
tribution describe the same number of particles, the temperature in the Tsallis 
45 distribution must be taken to be somewhat below that of the thermal distribu¬ 
tion. We will show this explicitly below. We will see that this change does not 
affect the conclusion that the rate of photon emission goes up and the times of 
emission are lengthened for a Tsallis distribution relative to the thermal one. 

In our analysis, we assume that the generalized (Tsallis) momentum distri- 
50 bution depends on only one variable, the temperature, i.e., the same principle 
scaling property as for thermal distribution functions. Such scaling behavior 
of distributions has been observed in simulations of the Glasma using classical 
field techniques [23 |2111301 [HI IHl 133] . It is of course more complicated if the 
longitudinal momentum and transverse momentum scales evolve differently in 
55 time. For simplicity, we assume only one scale, but the analysis could be ex¬ 
tended. It would of course be best to take the distribution directly derived from 
a first principle computation, and this might be possible in the future. 

Note that the distribution Q tends to 1 as i? —?► 0. At high px on the other 
hand, this distribution scales as {T/E)°‘. The dependence of the multiplicity 
upon the number of participants at any time t is 

where A is the number of participating nucleons. So for low energies, where 
/ —^ 1, the multiplicity distribution scales as the number of participants. 

60 However, at high energies, it scales as In saturation models, the 

initial temperature scales like T ^ A^/^ so that for a ~ 6, we get a very rapid 
^5/3 growth in the multiplicity0 

Note also that the low momentum part of the distribution does not evolve 
very rapidly in time while the high momentum piece falls as so that for 

^To understand this dependence, note that the density of partons in the transverse plane 
p determines the saturation scale p 
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3 REPLACING THERMAL DISTRIBUTIONS BY TSALLIS 

DISTRIBUTIONS 


&5 a = 


6 it fall as / ~ 1/t^ at fixed E. This is however less rapid than e 


for the Boltzmann distribution. This means there can be more radiation at later 
time. 

These considerations suggest that introducing power law tails for quark and 
gluon distributions can enhance the photon radiation rate and allow the radia- 
70 tion to appear at later times. At later times more flow will have been built up 
and the produced photon spectra will reflect that. Hence, this mechanism has 
the potential to solve the photon flow problem discussed in the introduction. 

To understand how the tails might affect the observed distributions of pho¬ 
tons, let us consider the radiation from an exponential distribution and compare 
75 it to that of a Tsallis distribution. For the following estimate we simply replace 
the photon distribution by a Tsallis distribution. We will improve on that by re¬ 
placing quark and gluon distributions by Tsallis distributions and recomputing 
the photon rate in the following section. 

Let us take the formula for thermal radiation Eq. Q ignoring logarithms 
and constant factors and integrate it over time: 



( 10 ) 


This assumes that the temperature of emission T ~ Ej/A is within the range of 
integration over temperatures. Now take a Tsallis distribution 



( 11 ) 


The stationary phase point of the integral is at 



( 12 ) 


For 0 = 6, this is 12, corresponding to a large change in the temperature of 
80 emission, which would make for a huge shift in the emission time, which goes as 
the cube of the temperature. Clearly such a big shift would move the emission 
outside of the range of integration over temperature where the QGP assumption 
is motivated, which will reduce this effect. Notice also that the value of the 
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4 DETAILED CALCULATION 


integrand at the stationary phase point is 256/9 which is much greater than the 
85 corresponding numerical factor for a Boltzmann distribution. 

This is a dramatic example of the effect of tails of distributions. Our more 
quantitative results in the following section are fortunately not so dramatic 
as this example. When we use Tsallis distributions for quarks and gluons to 
compute the rate of emission of photons by scattering, the result is indeed 
90 flatter than when using the thermal distributions. However, the enhancement 
in the kinematic region of emission for photons is not as large as in the previous 
example. Nevertheless we find significant enhancements of the photon yield and 
large lengthening of the photon emission times. Another factor that will reduce 
the effect is the overall lowering of the temperature for a Tsallis distribution 
95 relative to a thermal one when keeping the multiplicity fixed. 

We also need to discuss how the Tsallis distribution for quarks and gluons 
can arise. In p+p interactions, it is thought that Tsallis like tails arise from 
hadronization of jets. In the QGP the particles we are describing have not had 
time to turn into hadrons. In fact they must be generated from particles with 
100 momenta of the order or below the saturation scale. A Tsallis distribution is a 
reasonable guess for the distributions in this region, because the initial produc¬ 
tion of mini jets naturally leads to power law tails in addition to softer thermal 
particles. Note that in this picture, the initial typical momentum scales will be 
proportional to the saturation momentum, and the initial time to ~ l/Qsi so 
105 that at least for some time in the evolution we expect that the distributions will 
scale with the saturation momentum. If this is the case, then the A^part depen¬ 
dence of the initial distribution in the tail of the Tsallis distribution will rise 
much more rapidly than the Ncoii assumed for hard particle scattering. Since 
the Hnal state distribution of hadrons produced in this kinematic region is not 
no dramatically enhanced, final state interactions, that is jet quenching, must be 
very important. 
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4 DETAILED CALCULATION 


4. Detailed calculation 


In this section we present a detailed numerical calculation of photon produc¬ 
tion from a quark gluon plasma, using general distribution functions for quarks 
115 and gluons. We compare the numerical results in the thermal limit to the ana¬ 
lytic solution in [11117], and present the effect of using Tsallis distributions on 
the photon rate and typical emission time. 


4-1. Photon production rate 

The emission rate of photons with four-momentum Q = {E^, q) is given by 


E. 


dR 




(13) 


from the trace of the (12)-element ni 2 = of the photon-polarization tensor. 
In a thermal system the rate takes on the form miiiisTiiiTj 


dR _ 2 T 1 

( 2^)3 


(14) 


with the retarded photon self energy 11^,^. It is valid to all orders in the strong 
120 coupling as and to leading order in aem- 

When approximating the photon self energy by carrying out a loop expansion 


to some finite order Eq. (131 becomes equivalent to a description in terms of 


relativistic kinetic theory. Generally, expanding the self energy up to L loops is 
equivalent to computing the contribution from all reactions of m particles going 
125 to n particles, with m n < L -\- 1, with each amplitude calculated to order 
g^~^. Cutting the one-loop diagram for the photon self energy gives zero for 
an on-shell photon since the process qq ^ j has no phase space. Certain cuts 
of the two-loop diagrams give order corrections to this nonexistent reaction. 
Other cuts correspond to the reactions qq —)■ 57 , qg —>■ qj and qg —>■ 97 , the 
130 annihilation and Compton scattering processes. 

The contributions of these processes to the photon production rate are 
dRj 




/l(kl) / /2(k2) / (l±/3(k3)) 


x{27tAS^'^\Ki +K 2 -K 3 - , 


(15) 
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4.2 Numerical evaluation 


4 DETAILED CALCULATION 


with = J d?ki/{2ki{2'KY). The fj are the appropriate distribution functions 
and there is either a Bose enhancement or Pauli blocking term depending on 
the nature of the strongly interacting particle in the final state. The pre-factor 
follows from degeneracy factors and a sum over the charges of u and d quarks 
135 and one obtains .^^Cjoinpton — 20 and A/Annihilation — 320/3. 

The matrix elements squared are 

=-327r^egaas^ -k 0 (16) 

for the Compton process and 

= 327r\laas^ (17) 


for the annihilation process. 


The matrix elements (16) and ([17|) have poles at t and/or u = 0, which causes 


infrared divergent production rates ( |15[ ). The screening of these divergences 
is taken care of by including many body effects in a hard loop resummation 
140 for soft momentum transfers p < p* as done in [261127j . where p = |p| and 
P = (w,p) = Ki — Q. In the thermal case and the limit g —>■ 0 the total rate 
becomes independent of the cutoff p* and one obtains Eq. Q. 


4-2. Numerical evaluation 

In this work we are interested in the effect of non-thermal quark and gluon 
145 distribution functions on photon production. Even in this case it may be possible 
to perform a systematic hard loop resummation as done for anisotropic plasmas 
in [35]. We leave this study to future work and concentrate here on the effect 
of the non-thermal distributions on the rates in the region p > p*. 

The general rates for the annihilation and Compton processes in this region 
150 read 




dR. 


dP’q 




/^(p + q) f /g(k) l-(-/g(p-|-k) 


|p + q| 


|p + k| 


X <5(w -I- fc - Ip -I- k|) Q{p-p*) 


. t 


(18) 













4.2 Numerical evaluation 


4 DETAILED CALCULATION 


and 


E^ 


dRc, 


tPq 


= -IQ-asa 


/^(p + q) f /g(k) 1 - /,(p + k) 


|p + q| 


|p + k| 


X S{uj + fc - Ip + k|) Q{p - p*) 


t 


(19) 


where we relabeled k 2 to k and combined the two terms in the matrix element 
for the annihilation process assuming that the distribution functions for quarks 
and anti-quarks are the same. 

We simplify these expressions following [5B] and evaluate them using Monte 
Carlo integration. 

To study the effect of non-thermal tails in the distribution functions, we 
compare the thermal rates obtained by using Bose-Einstein and Fermi-Dirac 
distributions for gluons and quarks, respectively, with rates where those distri¬ 
butions were replaced by their Tsallis counterparts: 

-1 

, ( 20 ) 


n{E,) = 


E.y^ 

a 

1+ ^ 

) +1 

aT y 

/ 

E.y'^ 

a 

1+ ^ 

) -1 

aT, 

/ 


( 21 ) 


where the parameter a characterizes the power law tail of the distribution. We 
show a comparison of the thermal rates to those defined in (20) and (HI) in 

Fig-B 

160 We present a comparison of the thermal photon rates with those obtained 
using Tsallis distributions with parameter a = 6 for quarks and gluons in Fig.^ 
We chose the infrared cutoff p* = 10“°'^® « 0.56 to lie in the region where the 
soft and hard contributions are approximately equal. This means the full rate 
(sum of the soft HTL resummed and the hard contribution) is approximately a 
165 factor of 2 larger. 

One can see that the effect of the tails in quark and gluon distributions is 
significant. The photon rate becomes larger and harder. To show the effect 
more clearly we plot the ratio of the Tsallis to the thermal rate for different 
parameters a in Fig.j^ In fact, the curves shown in Fig. B are exponential fits to 
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4.2 Numerical evaluation 


4 DETAILED CALCULATION 



100 


0 2 4 6 8 


4 

E/T 


Figure 1: Comparison of thermal and Tsallis distribution functions with a = 6. E/T is the 
energy of the quark or gluon scaled by the temperature T. 

170 the ratios. We find that using Tsallis distributions with parameter a = 6 leads 
to almost exponential photon spectra but with an effective temperature twice 
as large as the temperature scale used in the calculation. One can also see in 
Fig.i that the Tsallis result approaches the thermal result very slowly. Only 
for parameters a « 1000 is the difference in rates negligible. 

Fig.|^shows that the ratios of Tsallis to thermal rates are almost independent 
of the infrared cutoff p*. This motivates us to modify the full thermal rate 
obtained using HTL resummation by the same factor to compute the total 
photon yield. Explicitly, for a = 6 we find that the ratio is well approximated 


by 



= 2.06+ 0.238 exp(0.608E.^/r). (22) 


To compute the produced photons from a one-dimensionally expanding quark 
gluon plasma, we integrate the photon rate using a time dependent temperature 
T{t) = To{to/ty/^, and set Og = 1/3. The transverse size of the system is fixed 


to At = -nB? with R = 6.6 fm, the radius of a lead nucleus. In the thermal 


case, we choose Tq = 500 MeV, starting at time <o = I/Tq. This value for Tq 
is extracted from numerical fluid dynamic simulation results at LHC energies 


10 












4.3 Photon emission times 


4 DETAILED CAL CULATION 



Figure 2: Comparison of thermal and Tsallis (a = 6) rates for p* Pi 0.56. Also shown is the 
analytic result for the hard part of the thermal rate from ED, which agrees perfectly with our 
numerical result. 


that produce hadron spectra and flow harmonics in agreement with experimen- 
tal data [39]. In case of the Tsallis distribution we multiply the thermal rate 
Q by the factor r{E^jT{t)) given in (22), and re-scale the initial temperature 
to ensure the same total number of quarks and gluons as in the thermal case. 
Assuming Boltzmann distributions, this is achieved when re-scaling the initial 
temperature in the Tsallis case by a factor 

1/3 


Nria) = 


( r(a) 


\a3r(a — 3) 


(23) 


175 which is Nt ~ 0.65 for a = 6. In case of the Bose or Fermi distributions, we get 
slightly different results of Nt{6) = 0.67, and Nxid) = 0.64, respectively. For 
simplicity, we will use the factor obtained assuming Boltzmann distributions in 
this work. The obtained photon yields are shown in Fig.[^ 


4-.3. Photon emission times 


We now compute the typical emission times for photons of different energies 
according to 


/to”"** {dR{t)/dy prdpr) 
tdt {dR(t)/dy pTdpx) 


(24) 


II 











4.3 Photon emission times 


4 DETAILED CALCULATION 



Figure 3: Exponential fits to the ratio of the Tsallis to the thermal rate for p* 5^ 0.56 and 
different a. Only for very large a does the Tsallis rate approach the thermal limit. 


(t) [fm/c] 

th. (ro = 500MeV) 

Tsallis (To = 325 MeV) 

PT = 1 GeV 

3.9 

4.6 

PT = 2 GeV 

1.7 

3.5 

PT = 3 GeV 

0.9 

2.5 


Table 1: Mean photon emission time using thermal and Tsallis rates and tmax = lOfm/c. The 
harder tails of the Tsallis distribution lead to significantly longer emission times compared to 
the thermal case, even when adjusting the initial temperatures to ensure equal numbers of 
partons in the system. 


180 Results using tmax = lOfm/c are presented in Table[^ We find that the harder 
tails in the Tsallis distributions lead to significantly longer mean emission times 
for photons. This effect can potentially delay photon emission so much that 
photons will be produced from a medium that has most of its anisotropic flow 
built up, leading to an increased photon elliptic flow. Detailed calculations of 
185 this using state of the art 3+1 dimensional hydrodynamics is the subject of 
future work. 

Finally, we show that the slope of the photon spectrum is quite sensitive 
to the time dependence of the temperature evolution. Decreasing the power 
1/3 by 10%, such that T{t) = To{to/t)^^^^, leads to a steepening of the photon 
190 spectra as shown in Fig.[^ In particular, this leads to better agreement with 
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5 SUMMARY AND CONCLUSIONS 



Figure 4: Dependence of the ratio of Tsallis to thermal rate on the infrared cutoff p* for a = 6. 


the experimental data in the low momentum region. Note that here we have 
not fixed the number of quarks and gluons to be the same as in the case of the 
usual 1+1 dimensional expansion, since we only want to emphasize the effect 
on the shape of the spectrum. 

195 5. Summary and Conclusions 

We have shown that a relatively small modification of thermal distributions 
to include power law tails results in dramatic modifications of the computed 
thermal photon yields and the times of emission. This modification has the 
correct properties to solve the problems with a theoretical description of photons 
200 from heavy in collisions as encountered for RHIC and LHC energies. To establish 
this explicitly, one must incorporate the modifications introduced in this work 
into a realistic 3+1 dimensional hydrodynamic simulation of the collision. 

To compute the power law tails from first principles may involve either a 
proper classical field computation, or cascade simulations where power law tails 
205 are dynamically generated. Such computations would allow one to assess the 
reasonability of our assumptions that the tails scale with the same momentum 
~ T, as does the center of the distribution, and that the power is close to a = 6. 
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Figure 5: Photon yields after emitting for lOfm/c and infinity in the thermal case with 
To = 500 MeV and using Tsallis quark and gluon distributions with Tq = 325 MeV and a = 6. 
The temperatures were chosen such that the number of quarks and gluons are the same in 
both cases. We compare to ALICE data for 0-40% central collisions from US]. 


The emission times we find are possibly long enough for part of the emission 
to occur in the hadron gas phase, where one may also expect power law tails 

210 of hadron distributions. It is however conceivable that in a strongly interacting 
medium such tails may be best described by partonic degrees of freedom as well. 
Such tails might have the property of geometric scaling since they arise from 
high energy processes, and thus would be consistent with the observed scaling 
behavior. 

215 Given the A dependence of the distribution functions in various transverse 
momentum regimes, our results suggest that the primeval distribution of quarks 
and gluons in heavy ion collisions has a significant Cronin enhancement at in¬ 
termediate transverse momentum values. This should have implications for jet 
quenching computations. 
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